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ABSTRACT 

We present an improved analysis of the final dataset from the QUaD experiment. Using an improved tech- 
nique to remove ground contamination, we double the effective sky area and hence increase the precision of our 
CMB power spectrum measurements by ^ 30% versus that previously reported. In addition, we have improved 
our modeling of the instrument beams and have reduced our absolute calibration uncertainty from 5% to 3.5% 
in temperature. The robustness of our results is confirmed through extensive jackknife tests and by way of the 
agreement we find between our two fully independent analysis pipelines. For the standard 6-parameter ACDM 
model, the addition of QUaD data marginally improves the constraints on a number of cosmological parameters 
over those obtained from the WMAP experiment alone. The impact of QUaD data is significantly greater for a 
model extended to include either a running in the scalar spectral index, or a possible tensor component, or both. 
Adding both the QUaD data and the results from the ACBAR experiment, the uncertainty in the spectral index 
running is reduced by ~ 25% compared to WMAP alone, while the upper limit on the tensor-to-scalar ratio is 
reduced from r < 0.48 to r < 0.33 (95% c.l). This is the strongest limit on tensors to date from the CMB alone. 
We also use our polarization measurements to place constraints on parity violating interactions to the surface 
of last scattering, constraining the energy scale of Lorentz violating interactions to < 1.5 x 10""*^ GeV (68% 
C.L). Finally, we place a robust upper limit on the strength of the lensing B-mode signal. Assuming a single flat 
band power between £ = 200 and £ = 2000, we constrain the amplitude of Z?-modes to be < 0.57 /iK^ (95% c.l.). 
Subject headings: CMB, anisotropy, polarization, cosmology 
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1. EVITRODUCTION 

Observations of the polarization of the cosmic microwave 
background (CMB) represent one of the most powerful probes 
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available for in vestigating the physics of the early universe 
(see e. g. Chall inor & Pei ris 2009 for a review). The CMB 
polarization field can be decomposed into two independent 
modes: even parity E-modes are generated, at the time of last 
scattering, by both scalar and tensor (gravitational wave) met- 
ric perturbations. In contrast, odd parity B-modes are gener- 
ated at last scattering only by gravitational waves, a generic 
prediction of inflation models. On small scales, Z?-modes are 
also expected to arise from gravitational lensing of the E- 
mode signal by intervening large-scale structures. A detection 
of B-mode polarization (on any scale) has yet to be made. 

After the initial detection of the much stronger /s-mode 
polarization (Kovac et ak 2002), steady improvements have 
been made in mea suring the £-mode signal by a num- 
ber of experiments ("Leitch et all 120051: iBarkats et al] 120051: 
Readhead et al. 2004; Montr ov et ail 120061: ISievers et all 
2007; Paaeetal. 2007; Wu et al.' '20071: lAde et alJ 120081: 
Bischoffetal. 2008; Nolta et al. 2009). Recently, a major 
step forward in precision CMB polarization measurements 
was achieved with the high-significance detection of a char- 
acteristic series of acoustic peaks in the E-mode polarization 
power spectrum with our initial analysis of the final QUaD 
dataset (Pryke et al. 2009; hereafter Paper II). For our analy- 
sis presented in Paper II, in order to mitigate against a strong 
polarized ground contaminant, we employed the technique of 
lead-trail differencing. Although this technique is extremely 
successful, it does have one major disadvantage — the effec- 
tive sky area is halved (while the signal-to-noise is kept the 
same) resulting in a corresponding increase of ^ 40% in the 
uncertainties on the final power spectrum estimates. The ma- 
jor improvement which we implement in this new analysis is a 
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technique to remove the ground contamination while preserv- 
ing the full sky area. Our analysis thus yields constraints on 
all six possible CMB power spectra which are approximately 
30% stronger than those presented in Paper II. 

We have also refined our modeling of the QUaD beams. 
For our previous analysis, we modeled the beams as ellipti- 
cal Gaussian functions. In addition to the main lobe, there 
is a small sidelobe component — with our increased sensi- 
tivity, we now find it necessary to explicitly model this side- 
lobe component. Accounting for the sidelobes results in a 
small 10%) increase in the amplitude of our power spec- 
trum measurements on small scales (multipoles, t > 700). 

Note that in this paper, we present the results obtained from 
two independent analysis pipelines. These are arbitrarily de- 
noted Pipeline 1 and Pipeline 2 and are derived from the two 
pipelines used to analyze the data from our first year of obser- 
vations (lAde et ai.ll20 08). The analysis presented in Paper II 
was performed using Pipeline 2. 

The paper is organized as follows. In Section|2] we briefly 
summarize the QUaD observations and low-level processing 
which are unchanged for this analysis. In Section [3] we 
present our new technique for mitigating against contaminat- 
ing ground pick-up and describe the details of our map mak- 
ing procedure. A description of our improved beam modeling 
is given in Section |4] and our treatment of the uncertainties 
is given in Appendix |A] The absolute calibration of QUaD 
is briefly described in Section |5] with an error analysis given 
in Appendix iBl The results from our two independent power 
spectrum analyses are presented in Section |6] In Section |7] 
we combine the QUaD results with data from the WMAP, 
ACBAR and SDSS experiments to place constraints on the 
parameters of a number of cosmological models. Our conclu- 
sions are presented in Section|8] 

2. OBSERVATIONS AND LOW-LEVEL PROCESSING 

The QUaD experimen t and its performance are described 
in lHinderks et al.l (l2009h . hereafter referred to as Paper I. The 
low-level data processing is described in Paper II. The initial 
low-level processing of the raw data has not changed for the 
analysis presented in this paper so here we give only a brief 
summary and refer the reader to Paper II for a detailed de- 
scription. 

The QUaD experiment was a 2.6m Cassegrain radio tele- 
scope which observed from the South Pole for three seasons 
from 2005 to 2007. The QUaD receiver consisted of 31 pairs 
of (orthogonal) polarization sensitive bolometers (PSBs), 12 
at 100 GHz and 19 at 150 GHz. These PSB pairs were ar- 
ranged on the focal plane in two orientation angle groups sep- 
arated by 45°. The raw time-ordered data (TOD), which were 
sampled at 100 Hz, were first deconvolved to correct for the fi- 
nite response times of the PSBs and electronics used to detect 
the incoming signal. The time-constants used for each detec- 
tor were measured using an external Gunn oscillator source 
as described in Paper I. After deconvolution, the detector data 
were low -pass filtered to < 5 Hz.' The data were then de- 
glitched to remove cosmic rays and other events. A relative 
calibration was then applied to each detector using the "eleva- 
tion nod" technique described in Paper I and Paper II. 

' QUaD's CMB observations employed a relatively slow scan speed of 
0.25°/ sec. For our observing declination (~ -50°), the sky signal for multi- 
poles, £ < 2000, appears in the time-stream at < I Hz. The low-pass filtering 
therefore removes none of the sky signal (for t < 2000) but it does remove 
high frequency noise introduced by the deconvolution procedure. 



Once de-glitched and calibrated, the data were downsam- 
pled to 20 Hz. For the analysis presented in this paper, we 
have retained the exact same data cuts for bad weather, moon 
contamination and badly behaved detectors as were used in 
Paper II. Out of a total of 289 days of observations during 
2006 and 2007, after applying these data-cuts, 143 remained 
for the science analysis. Although fully code independent, 
the low level parts of our two analysis pipelines are algorith- 
mically similar 

3. MAP-MAKING USING GROUND TEMPLATE REMOVAL 

Our improved technique for removing the ground signal re- 
lies on redundancies in the scan strategy so before describ- 
ing our templating procedure, we turn first to the QUaD scan 
strategy and examining the redundancies present within it. 

3.1. QUaD observing strategy 

QUaD observed a ^ 100 square degree area of sky, centered 
on RA 5.5h, Dec -50°. The field is fully contained within the 
shallow field observed by the 2003 flight of the Boomeran g 
experiment (hereafter referred to as B03, iMasi et aTl 120061) . 
The QUaD field also partially overlaps with B03's deep 
field. The QUaD observations employed a lead-trail scheme, 
whereby each hour of observations were split equally between 
two adjoining subfields, separated in RA by 0.5 h — the lead 
field, centered on RA 5.25 h, and the trail field, centered on 
RA5.75 h. 

The scanning strategy consisted of constant-elevation scans 
back and forth over a 7.5 deg throw in azimuth, applied as 
a modulation on top of sidereal tracking of the field center. 
Each hour of observation was equally split between the lead 
and trail fields. These half-hour sessions were further divided 
into four "scan-sets", consisting of ten "half-scans" each, and 
the telescope was stepped in elevation by 0.02 degs between 
scan-sets. After a half hour scanning the lead field, the tele- 
scope pointing returned to its starting position in azimuth and 
elevation, and repeated the same scan pattern with respect to 
the ground, but now scanning the trail CMB field. The trail 
field's scan pattern was thus a replica (in azimuth/elevation 
coordinates) of the lead field's. After an hour the pointing 
moved on to a fresh part of sky and the process repeated. This 
scan pattern was designed to facilitate the lead-trail differenc- 
ing analysis presented in Paper II, whereby each pair of lead- 
trail partner scans are point by point differenced. Any ground 
signal, which is stable in time over the half hour which sep- 
arates the lead and trail observations, will be completely re- 
moved by this differencing, at the expense of a reduction in 
the effective sky area by a factor of two.^ 

In addition to the lead-trail scheme, a further redundancy is 
present in the scan-strategy due to the movement of the CMB 
field across the sky during each scan-set. The time elapsed 
from when a given sky pixel is first visited on the first half 
scan of a set, to when it is last visited on the tenth half-scan 
is ~ 6 minutes. During this time the sky rotates by 1.5 de- 
grees. Using only data from a single scan-set, there is there- 
fore scope for separating signals originating on the ground 
from those originating on the sky, on scales smaller than 1 .5 

^ For a Gaussian field, such as the CMB, neglecting correlations between 
the signal in the lead and trail fields, the fluctuations in the differenced field 
will be amplified by a factor of \/2 and the power spectrum will increase by 
a factor of 2. The noise is amplified in a similar manner and so the signal-to- 
noise ratio in the differenced field remains unchanged from that achieved in 
the non-differenced field. 
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degrees, corresponding to ^ 250 in multipole space. One 
can achieve further separation of ground from sky by com- 
bining the data from lead and trail partner observations. For 
this work we have made use of both the lead-trail and sky 
rotation redundancies to separate the ground and sky signals, 
and to reconstruct the CMB fields over the full sky area. 

3.2. Ground template removal 

Field differencing is a sub-optimal use of the redundancy in 
the scan strategy to mitigate against ground pickup. One can 
retain more of the sky information by constructing and remov- 
ing estimates of the ground signal. To facilitate the removal 
of ground signal, we can model the TOD as 

^^,■ = 5,(61) + ^(a) + n,+ Oscan, (1) 
where 5,(0) is the sky signal, 

Sm = ^ [/(0) + e(0)cos(2(/.,) + f/(0)sin(20,)] . (2) 

Here, /, Q and U are the Stokes parameters in the direction, 6 
on the sky and 0, is the polarization sensitivity angle (a combi- 
nation of detector orientation on the focal plane and boresight 
rotation) of each detector 

In equation ([T]i, g{a) represents the ground signal as a func- 
tion of azimuth, a. In what follows, we will be constructing 
and removing estimates of g{a) for each QUaD detector and 
for each pair of lead and trail scan-sets independently. Within 
these subsets of the data, the elevation is constant and so the 
ground signal estimates (which we refer to as ground "tem- 
plates") are constructed as a function of azimuth only. How- 
ever, since we allow the templates to differ between detectors 
and between lead-trail scan-set pairs, in practice, the ground 
signal which we remove from the data does depend on eleva- 
tion also. Moreover, since templates are constructed for each 
detector individually, our model also allows the ground signal 
to depend on frequency. 

The resolution with which to construct the ground tem- 
plates, in general, needs to be determined through trial and 
error. In practice, we find that the effect of changes in this 
parameter on the resulting maps are imperceptible. This sug- 
gests that the ground signal varies smoothly in azimuth and 
that a fairly coarse resolution is sufficient to characterize it. 
For this analysis, we have used a resolution of Aa = 0.1 degs, 
midway between the minimum and maximum resolutions we 
have investigated. 

In equation ([1]), we have split the noise component into a 
random part (n,-, which we model as a Gaussian random vari- 
able; see Section l672l i and an offset (Oscan) which we model as 
a constant for each half-scan. We have found that it is essen- 
tial to remove these offsets before constructing ground tem- 
plates from the data. Otherwise the resulting templates tend 
to be dominated by the long-timescale part of the 1/f atmo- 
spheric noise (i.e. the offsets) rather than the ground signal 
which we are attempting to characterize. 

Regardless of the technique employed to reconstruct maps 
of the Stokes parameters (naive, maximum-likelihood etc.), 
one can re-cast the well known map-making equation (e.g. 
iStompor et alj|20()2l and references therein). 



by making the substitutions. 



m 

A 



m+g 

■(A^^A^), 



(4) 



where m is the reconstructed CMB map, g is the reconstructed 
ground signal and A'^'^* and A'^ are the "pointing matrices" 
associated with the CMB and ground signals respectively. For 
a highly redundant scan strategy, equation (O should be solu- 
ble exactly and the CMB and ground signals should be com- 
pletely separable. However, for a scan strategy such as that 
used for QUaD with limited revisiting of the same sky pixels 
at different azimuths, this complete separation between sky 
and ground is not possible. In this case, one can still separate 
the ground and sky signals on smaller scales but one loses all 
information on the largest scales where the separation is de- 
generate. For QUaD, we find that below £ ^ 200 our tem- 
plate removal procedure offers essentially no improvement 
over field-differencing and that the separation is near perfect 
by 1000. 

Attempting to simultaneously solve for the sky and ground 
signals in the QUaD data by applying equation (O, we find 
that large scale (ground signal) gradients are introduced to 
the resulting CMB maps due to the degeneracy between the 
CMB and the ground signal on the largest scales. One could 
certainly modify the procedure (e.g. by marginalizing over 
the large-scale CMB and ground-signal modes) to solve this 
problem. For the analysis presented here, we adopt a sim- 
pler approach and simply solve for the ground signal inde- 
pendently, subtract this from the TOD, and then construct the 
CMB maps from the ground-cleaned TOD. We account for 
the resulting filtering of the CMB signal in our Monte-Carlo 
analysis. Our analysis assumes that the ground signal does 
not change between the start of the lead scan set and the end 
of the corresponding trail scan set 36 mins). This is only a 
slight relaxation of the assumption that was made for our pre- 
vious analysis where we assumed that the ground signal was 
constant over a 30 minute timescale. 

To apply the template removal, we proceed as follows. 
First, we estimate and remove the atmospheric 1 // offsets, 
Oscan from each half scan. To estimate the offsets, we simply 
take the mean of the data within each scan. Note however 
that we restrict the azimuth range over which we calculate the 
offsets to the central azimuth range where all the scans within 
a scan-set overlap. This ensures that our estimated ground 
templates will be unbiased over the full azimuth range. 

After removing the offsets, templates of the ground signal 
are constructed by simple binning of the timestream data in 
azimuth. That is, for each ground template "pixel", we con- 
struct 

^("^=^E'^'' ^5) 

where the sum is over all data from the lead scan set and its 
partner trail scan set which falls in the azimuth range, Aa. 

Although they will be unbiased estimates of theground sig- 
nal, the templates constructed using equation (|5| will also 
contain CMB signal and noise. The expectation value of the 
constructed templates is 



A^N"'Am = A^N"'d, 



(3) 



-'Vhits 



J2 Si(0)+ J2 



to reconstruct both the sky signal and the ground signal simply 



(6) 

JeAa ieAa . 

where Si(9),g(a) and «,■ are the quantities defined in equa- 
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Fig. 1 . — Demonstration of the performance of the ground templating procedure described in the text. The figure shows maps of the Stolces U polarization over 
the full QUaD sky area at 150 GHz, smoothed with a 5 arcmin Gaussian kernel. The map on the left is that obtained without removing ground templates and 
is heavily contaminated by ground pickup. Note the similarity of the contamination between the lead and trail halves of the map. The map on the right is that 
obtained when we include our templating procedure. Clearly the vast majority of the ground signal is successfully removed by this process. For the purposes of 
this illustration, in order to highlight the success of the template removal, in both cases, only the mean of each scan was removed from the TOD before mapping. 
However, for our cosmological analysis, we use maps which have had third-order polynomials removed from each scan (see Figure|2j. 



tion ([T]i. In the limit of a highly redundant scan strategy and 
uncorrelated noise, the terms in brackets will average to zero 
and the templates will contain only ground signal.^ For a re- 
alistic scan strategy and correlated noise, these terms will be 
non-zero and so removing the templates will have the effect 
of filtering both large-scale CMB signal and long-timescale 
noise. We correct for this filtering by including the tem- 
plate removal procedure in our Monte-Carlo simulations (Sec- 
tion |62]i. 

Finally, to obtain estimates of the TOD which are free of 
ground pickup, the templates are subtracted from the original 
TOD; 

d'=''="" = d-A^g, (7) 

where A*^ is the pointing matrix associated with the ground 
signal and g is the estimated ground signal constructed using 
equation (|5]l. This will result in TOD which is, in principle, 
free of ground contamination and can thus be modeled as in 
equation ([TJ but now without the g{a) term. Once the ground 
contamination has been removed, our map-making proceeds 
as described in Paper II. Explicitly, we perform the follow- 
ing operations. For each azimuth scan the best-fit third order 
polynomial is subtracted to remove the long timescale part of 
the atmospheric 1 // noise. For each PSB pair, the data is then 
summed and differenced to yield temperature (5,) and polar- 
ization {di) TOD. Maps of the temperature (or Stokes /) CMB 
field are then constructed from the summed data using a sim- 
ple weighted average: 
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(8) 



where the weights are given by w; = W(jic)/vscan- Here, Vjcan is 
the variance of the data across the parent half-scan — noisy 
data (e.g. due to bad weather) is thus down-weighted. W(x), 

^ For an experiment sensitive to absolute temperature, the first term in 
brackets in equation (6) would, in fact, average to the CMB monopole for a 
highly redundant scan strategy. For an experiment sensitive to temperature 
differences only (such as QUaD), this term would average to zero. 



where x denotes the fractional position within the scan, is an 
apodization (the same for each scan) which we use to down- 
weight the scan ends. We apply this apodization to reduce 
the tiling effects seen in our previous analysis (see Section 
6.3 and Figure 15 in Paper II) whereby the interaction of the 
polynomial filtering with different sky coverage for different 
detectors produced visible step features in the final maps. The 
exact form used for the apodization is not important. 

We construct maps of the Stokes polarization parameters, 
Q and U as 



(cos^(2(/),-)) (cos(2(/);) sin(20,-)) 
cos(20,)sin(20,)) (sin^(20,)) 

^ cos(24>i)d, 
'sin(2(f)i)di 



(9) 



where the angled brackets denote an average taken over all 
data falling within each map pixel and the angle, 0, is a com- 
bination of the polarization sensitivity direction of each de- 
tector on the focal plane and the "deck angle" (rotation about 
the telescope boresight) of the observation. Note that to con- 
struct the averages (e.g. (cos(20,)t/,)) on the right-hand side 
of equation (|9]), we also use inverse-variance weights as in 
equation ([8]). 

Note that the only difference between the approaches of 
our two pipelines in applying the above operations is during 
the final coaddition of the template-subtracted data into CMB 
maps: Pipeline 1 uses HEALPix"* to pixelize the sky at a res- 
olution of ^ 1.7 arcmin (A^side = 2048). Pipeline 2 works un- 
der the flat sky approximation and pixelizes the sky into a 2D 
cartesian grid with a spacing of 1.2 arcmin. Figure [T] shows 
an example of the performance of the template removal pro- 
cedure (for Pipeline 1) for the 150 GHz Stokes U polarization 
map. Note that, in order to highlight the success of the tem- 
plate removal, for this demonstration we have not applied the 
third order polynomial removal mentioned above to the TOD 

See http://healpix.jpl.nasa.gov/index.shtml| and lG6rski et al] 42003) . 
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and have only removed the mean from each scan.^ 

In Figure |2] we present the full set of maps {T,Q and U 
at 100 and 150 GHz; again for Pipeline 1) over the full sky 
area as estimated using the template removal procedure (now 
including the third order polynomial removal). For the pur- 
poses of visual illustration only, we have smoothed each of 
the maps with a 5 arcmin Gaussian kernel in order to bring 
out the CMB structure. 

We can quickly (and crudely) assess the relative amounts 
of E- and B-mode power in the polarization maps by decom- 
posing the Q and U maps into E and B modes. We do this 
under the flat-sky approximation. To minimize the impact of 
the noisy edge regions of the maps, and to reduce the effects 
of E/B mixing due to the finite survey geometry, we apply an 
apodization to our maps before Fourier transforming. (E/B 
mixing is fully accounted for during our power spectrum es- 
timation described in Section|6]) The resulting E and B maps 
at 150 GHz, are shown in Figure [3] We clearly detect sig- 
nificantly more E-mode than Z?-mode structure. The recon- 
structed B-mode map shows similar levels o f flu ctuations to 
our polarization jackknife maps (see Section |63] | and is con- 
sistent with noise. 

4. BEAM MEASUREMENTS AND MODELING 

Our main set of observations for investigating the QUaD 
beam shapes are a series of single day observations of the 
QSO PKS0537-441. In Paper II, our beam model consisted 
of elliptical Gaussian fits to these QSO observations for each 
channel. 

Given our increased sensitivity, we now include an addi- 
tional sidelobe component in our beam model. In order to 
measure the sidelobes from the QSO data, we apply a 6th- 
order polynomial filter to the TOD before mapping (with the 
QSO masked) and coadd these data over all channels and over 
three days of observations. The radially averaged beam pro- 
files measured from these maps reveal the presence of side- 
lobe structure at just below the -20 dB level, as predicted by 
the physical opt i cs (PO ) simulations of QUaD presented in 
lO' Sullivan et al.l (12008 1). Our two analysis pipelines model 
these observations in slightly different ways though both are 
matched to the QSO data. 

Pipeline 1 rescales the PO models. The beam profiles which 
directly result from the PO simulations are not a perfect match 
to the QSO observations. In particular, the predicted main 
lobe widths are smaller than observed while the predicted lev- 
els of sidelobes are somewhat larger than observed. To match 
the PO models to the QSO data, we parametrize the models 
using two parameters, one which scales the main lobe width 
and one which varies the amplitude of the sidelobes. We then 
use the observed QSO radial profiles to fit for these param- 
eters. The resulting best-fit re-scaled PO models are used to 
model the beam. 

Pipeline 2 models the beams in a fully empirical manner 
and is an extension of the model used in Paper II. Using 
the existing elliptical Gaussian fits to the quasar data, a pure 
Gaussian simulated beam coadded across channels and obser- 
vation dates is generated and subtracted from the measured 
QSO maps. The residual after subtraction is the sidelobe com- 
ponent of the beams. This residual is too noisy to be used di- 
rectly and so it is radially averaged to produce an (assumed) 

^ In fact, the polynomial fitting procedure does remove the gross features 
of the ground signal from the data although much remains — full field maps 
which have not been subjected to template removal fail jackknife tests at high 
significance regardless of whether we apply polynomial removal or not. 



azimuthally symmetric sidelobe template. This sidelobe tem- 
plate is then added to the original Gaussian elliptical models 
to produce a fully empirical beam model. 

Figure |4] shows the radially averaged profiles measured 
from the QSO data along with the profiles as predicted us- 
ing our old elliptical Gaussian beam model and as predicted 
using our current beam models. Our revised beam models 
are clearly a superior description of the true beams and are 
in good agreement — in terms of the resulting beam trans- 
fer functions, the two beam models agree to within 4% at 
100 GHz and to within 2% at 150 GHz for £ < 2000. A de- 
scription of how we account for the remaining uncertainties 
on our beams is given in Appendix lAl 

5. ABSOLUTE CALIBRATION 

We derive the absolute calibration of QUaD by cross- 
correlating our temperature maps with maps from BOS 
(iMasi et al.l 12006 ). This analysis is done in spherical har- 
monic (aem) space following the calibration technique used by 
Boomerang wh ich, in turn, was calibrated against the WMAP 
Ist-year maps (iBennett et al. 2003).* We apply a correction 
to the B03 maps to account for the change in calibration 
(1.25% in temperature) between the 1- and 5-year WMAP 
analyses (Hinshaw et al. 2009). The B03 maps (which are 
essentially unfiltered) are first passed through the QUaD sim- 
ulation pipelines to ensure that they are filtered in an identical 
manner to the QUaD maps. Taking the spherical harmonic 
transforms of the maps, the absolute calibration factors for 
QUaD are then given by 

nquad / b03-X * b03-Y\ 
yquad _ "e \"em "em I qq\ 

'^l \"lm °-lm I 

where the superscripts, X and Y denote two noise-independent 
145 GHz B03 maps and Z?^"'"' and B^"^ are the beam transfer 
functions of the two experiments. This process produces an 
absolute calibration factor (in /iK/V) as a function of multi- 
pole, £. We take our final absolute calibration factors to be the 
mean of this value in the range 200 <l < 800, corresponding 
to the overlapping scale range of the two experiments. The 
largest contributors to our calibration error are the quoted B03 
uncertainty of 2% and a relative pointing uncertainty between 
QUaD and B03. Our final calibration uncertainty is 3.4%. 
Further details on how we estimate this uncertainty are given 
in Appendix iB] 

6. POWER SPECTRUM ANALYSIS 

To estimate the CMB power spectra, as in Paper II, we 
adopt a Monte-Carlo (MC) based technique whereby we rely 
on accurate simulations of the experiment to correct for the ef- 
fects of noise, beams, timestream filtering and the removal of 
the ground templates. Before describing the MC simulations, 
we first describe the differences between our two pipelines in 
their approach to power spectrum estimation. 

6.1. Power spectrum estimation 

Both of our pipeUnes broadly f ollow the so-called pseudo- 
Ct technique (|Hivon et al.l l2002h . extended to polarization 
dBrown et al.ll2005l) . Note that, for both pipelines, in order to 

* We have also performed the calibration against the WMAP maps directly 
and we find consistent results. However, calibrating using the B()3 maps 
produces a more accurate result due to the larger overlap in angular scale of 
B03 and QUaD. 
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Fig. 2. — Maps of the T (top panels), Q (middle panels) and U (lower panels) Stokes parameters over the full QUaD sky area at 100 GHz (left) and 150 GHz 
(right). For display purposes only, the maps have been smoothed with a 5 arcmin Gaussian kernel. Note the difference in the color stretch used to display the 
temperature and polarization maps. 
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Fig. 3. — Apodized 150 GHz QUaD polarization maps decomposed into £-modes (left) and B-modes (right). Once again, the maps have been smoothed with 
a 5 arcmin Gaussian kernel. The data is clearly dominated by fi-modes. The amplitude of fluctuations present in the B-mode map is consistent with noise. The 
sUght reduction in the amplitude of fluctuations towards the central RA of the field is due to the application of the apodization mask which down-weights the 
"seam" between the lead and trail halves of the map. 




Fig. 4.— QUaD beam profiles at 100 GHz (top panel) and 150 GHz (bot- 
tom panel) as measured from the QSO PKS0537-441. The radial profiles as 
predicted using our new beam models for both pipelines are over-plotted as 
the red and blue curves and show good agreement with the QSO data. Also 
shown for comparison is the elliptical Gaussian beam model used for our 
previous analysis in Paper 11. 

minimize edge effects, the T, Q and U maps are first apodized 
with an inverse-variance mask as described in Section 6.1 of 
Paper II. 

Pipeline 1 works on the curved sky and uses fast spherical 
harmonic transforms to estimate the pseudo-Q spectra. These 
spectra are then corrected for the effects of the sky cut, noise 
and filtering, and binned into band powers according to 

Pi = ^K^^, ^P,-,(Q-(N,)mc). (11) 



In this equation and throughout this section, we use bold-face 
to denote six-spectrum quantities, e.g. is a 6 n band vector, 
Pb = {Pl^Pl^^Pr^PrA^^Pf,"}^ and Kw is a 6«band x 
6 n band matrix, given by 



Kfei- = Pbe J2 Qi'b' 



(12) 



In equation (fTTl l. are the pseudo-Q spectra estimated from 

the data maps and (N^ )mc are the noise power spectra as mea- 
sured from simulations. 

Phi is a binning operator which bins the raw pseudo-Q into 
band powers and Qa, is the inverse operator which "unfolds" 
a band power into individual Qs. For this analysis, we use 
"flat" band powers for which the quantity, £{£+l)Ce/2TT is 
constant within each band. That is, the binning operator we 
use is 



Phe = 



■^low 

0, 



if2<C<^< 
otherwise. 



«(i'+i) 

"-low 



with the corresponding inverse operator given by 

if2<4'i<^<Aow^ 



Qt 



0, 



£+1) ' 



low 

Otherwise, 



(13) 



(14) 



where if^^ denotes the nominal lower edge of band b. 

The coupling matrix, M^^/, describes the mode-mixing ef- 
fects of the apodiz ation mask and sky cut and is given in 
iBrown et al.l (12005b for the full set of six possible CMB spec- 
tra. Note that the matrix fully encodes E /B leakage ef- 
fects due to the finite survey geometry and so our Pipeline 1 
estimator explicitly corrects the EE and BB spectra for this 
leakage in the mean. 

/v is a transfer function which we use to describe the com- 
bined effects of timestream filtering, beam suppression and 
filtering of the sky signal due to the removal of the ground 
templates. This function will also encode any other signal 
suppression effects which are present in our simulations (e.g. 
pixelization effects). We estimate from our signal-only 
simulations as described in Section l63] 

Pipeline 2 works in the flat sky approximation and uses 2D 
FFTs to estimate power spectra. This pipeline is described in 
detail in Paper II. The power spectrum estimator for Pipeline 
2 can effectively be written as 



(15) 
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where Fb is the binned equivalent of the per-multipole trans- 
fer function, F( and we have implicitly made the connection 
between the flat sky and curved sky power spectra, Q w P{k) 
for I « 2-Kk. 

Note that (in addition to the flat sky approximation), the pri- 
mary difference between the two pipelines is that Pipeline 1 
performs the correction for the mode-mixing effects induced 
by the sky cut whereas Pipeline 2 does not perform this cor- 
rection. Because of this difference, neither the recovered 
band powers nor their uncertainties are directly comparable 
between the two pipelines. A proper comparison of the two 
analyses requires the u se of the associated band power win- 
dow functions (Section 16.3b which fully encode the relation 
between underlying true sky power and observed power for 
both pipelines. 

In both analyses, we estimate the covariance matrix of our 
power spectrum estimates from the scatter found in the power 
spectra measured from simulations containing both signal and 
noise: _ _ 

(APfcAPiO = {{Vb-Vb){Vb'-n'))Mc. (16) 

where denotes the average of each band power over all 
simulations. Note finally that the covariance properties of the 
power spectra estimates are dependent on whether the correc- 
tion for mode-mixing induced by t he s ky cut is applied or not. 
We return to this issue in Section l676l where we compare the 
results from our two analyses. 

6.2. Simulations 

In simulating QUaD, we follow the procedure described 
in Section 5 of Paper II with some important differences, 
which we now discuss. Algorithmically both of our analy- 
sis pipelines adopt the same approach to creating simulated 
timestreams and onl y dif fer in the final map-making stage as 
described in Section [J!2] 

6.2.1. Signal simulations 

To create the signal component in the simulations, we first 
generate m odel TT,EE,TE a nd BB CMB power spectra us- 
ing CAMB ( iLewis et alj2000h . The input cosmology consists 
of the best-fitting ACDM model to the 5-year WMAP data set 
(|Dunkley et al. 2009). Note that the model spectra used in- 
clude the effects of CMB lensing and so the input B-mode 
power is non-zero. (For comparison, in Paper II, our input 
model was the best-fit to the 3-year WMAP data set and our 
input B-mode power was set to zero.) 

Realizations of CMB skies are then generated from these 
model spectra using a modified version of the HEALPix soft- 
ware. These maps are generated at a resolution of 0.4 arcmin 
(A^side = 8192). The simulated maps are then projected onto 
a 2D cartesian grid and convolved with the beam model for 
each detector channel. The resolution used for this interme- 
diate map is 0.6 arcmin. The generation of the sky maps for 
each detector and deck angle and interpolation to simulated 
TOD then proceeds exactly as described in Section 5. 1 of Pa- 
per II. 

6.2.2. Noise simulations 

In order to simulate realistic noise, we must first measure 
the noise properties from the real data. However, the undiffer- 
enced data contains not only noise but also CMB signal and 
ground signal (see equation [T]i. The instantaneous signal-to- 
noise in the timestream is negligible and so the CMB com- 
ponent can be safely ignored. However, the same is not true 



for the ground signal which, in some cases, is a very signif- 
icant component in the timestream, particularly in polariza- 
tion. The data must therefore be cleaned of the ground com- 
ponent before measuring the noise power spectra. 

One might think that the best way to achieve this would be 
to measure the noise spectra from data which has been tem- 
plate subtracted using equation (|7]l. Such a procedure could, 
in principle, be iterated and is similar to procedures suggested 
for measuring the noise properties of CMB data when the sig - 
nal component is non-neglig ible (e.g.'Ferreir a & JaffdllOOOl) . 
However, as noted in Section [J!2l our template removal proce- 
dure filters the noise in a non-trivial fashion. In particular, be- 
cause of the non-uniform azimuth coverage of the scan strat- 
egy, the ground templates are noisier at each end than they are 
in the central regions. The result is that after subtracting the 
templates, the noise is no longer uniform and this prohibits its 
characterization through simple FFT-based power spectrum 
estimators — since the noise is no longer a stationary Gaus- 
sian random process, a power spectrum description will fail. 

To avoid these complications, we have measured the noise 
properties of the TOD from data which has been lead-trail 
differenced. The differencing efficiently removes the ground 
signal while under the assumption that the noise is stationary 
over a 30 minute timescale, the power spectra of the undiffer- 
enced TOD are simply the spectra measured from the differ- 
enced TOD divided by 2. For each pair of lead-trail observa- 
tions, we therefore assign the power spectra measured from 
the differenced data, appropriately normalized, to each of the 
lead and trail scan-sets. Simulated noise-only timestreams are 
then generated exactly as described in Section 5.2 of Paper II. 

Examining the QUaD data and comparing it to simu- 
lated data obtained using the above process, there are occa- 
sions where our assumption of stationarity over a 30 minute 
timescale is not satisfied. However, for the majority of the 
data the assumption is good and it is only ever a poor one for 
our temperature analysis. Moreover, a thorough comparison 
of the statistics of the simulated and real data indicates that 
our procedure provides an excellent description of the noise 
properties of the undifferenced data when averaged over each 
day for both temperature and polarization. Further averag- 
ing over tens of pixels and hundreds of observation dates will 
result in these rare failures of our noise model having a negli- 
gible impact on the results. 

Once generated, both the signal-only and noise-only simu- 
lated TOD are processed into T, Q and U maps in an identical 
manner to that used for the real data. In particular, note that 
both the ground template subtraction and polynomial filtering 
are applied also to the simulated data and so the effects of fil- 
tering on both the signal and on the noise are fully accounted. 
Finally, to obtain simulated maps containing both signal and 
noise, we simply add the signal-only and noise-only maps. 
Since all of our data processing steps are linear operations, 
this final step results in simulated maps no different to those 
which would have been obtained if we had instead summed 
the signal-only and noise-only TODs, and is computationally 
more efficient. 

6.3. Transfer functions and band power window functions 

We estimate the transfer function from our suite of signal- 
only simulations. In the absence of noise, the mean of the 
recovered pseudo-Q spectra will equal their expectation val- 
ues, 

(Q)mc = 5]m«'F,,C,., (17) 
I' 
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Fig. 5. — Total transfer function, Fg (black curve) as measured in our 
Pipeline 1 analysis for the 150 GHz channel. Also shown are the transfer 
functions for timestream filtering (green), ground-template removal (blue) 
and beam suppression (red) in isolation. The dashed curve shows the suppres- 
sion of signal due to map pixelization for a HEALPix resolution, Wside = 2048. 
The total transfer function is the product of these four individual curves. 

where are the input model spectra used to create the simu- 
lations. For a small-area survey such as QUaD, the unbinned 
coupling matrix, M^^/, is singular and so equation ( [TtI i can- 
not be solved directly. In Pipeline 1, we iteratively solve this 
equation to provide an estimate of F^. With a reasonable start- 
ing guess, convergence is typically reached in just a few itera- 
tions. For Pipeline 1, the band power window functions, Wa£, 
defined by (.Knox 1999), 



27r 



are given by 

I 



1-K 



(18) 



(19) 



Pipeline 2 calculates its band power window functions nu- 
merically as described in Section 6.6 of Paper II. In order to 
calculate the transfer functions. Pipeline 2 simply takes the ra- 
tio of the mean band powers recovered from signal-only sim- 
ulations and the expectation values for each band power: 



(20) 



Figure |5] shows the derived transfer function from Pipeline 
1. (The transfer function from Pipeline 2 is similar) As 
mentioned earlier, this function encapsulates all effects due to 
timestream filtering, beam suppression and the filtering due to 
removal of the ground templates. To demonstrate the relative 
size of these effects, we also plot the transfer function derived 
from special simulations with these three effects included in 
isolation. Of particular interest is the transfer function de- 
scribing the ground-removal procedure — this curve, in ef- 
fect encapsulates the lossiness of the technique. On scales 
where this curve is > 0.5, we gain an improvement over 
the field-differencing technique. We see that below ^ ^ 200, 
there is little gain from our template-removal technique over 
field-differencing. The amount of signal retained then climbs 
rapidly and is effectively unity by ^ = 1000. 
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Fig. 6. — Band power window functions for Pipeline 1 (top panels) and 
Pipeline 2 (bottom panels) for the C^® power spectrum. The blue curves 
show the response to true BB power and the red lines show the response to 
EE power 

Figure |6] shows the band power window functions (BP- 
WFs), plotted as Whijt, for from both our pipelines for 
the 150 GHz channel. These functions describe the response 
of our band power measurements to the true sky signal at each 
multipole. The negative wings in the BPWFs of Pipeline 1 are 
a direct result of the application of the correction for mode- 
mixing effects as described in Section |67T] 

The expectation values for our EE band powers assuming 
the best-fit ACDM model to the WMAP 5-year data are shown 
in Figure [T] Note that the apparent improvement in going 
from Pipeline 2 to Pipeline 1 in terms of the agreement be- 
tween the true sky power and the band power expectation val- 
ues does not come without a price — as a result of applying 
the mode-mixing correction, the error bars for Pipeline 1 are 
enhanced with respect to the Pipeline 2 errors. The covariance 
properties also change between the two analyses such that the 
total information content is preserved. 

We emphasize that, in terms of either the accuracy or the 
precision of the recovery of true sky power, neither analysis 
is superior. This is clear from the fact that one can transform 
between the band power estimates and covariances of the two 
analyses via a simple and exact matrix operation. Whether to 
apply the mode-mixing correction is simply a matter of pref- 
erence and is only relevant for visual interpretation of the re- 
sults — one can choose to have smaller error bars or one can 
choose to have band powers which better trace the underlying 
true sky power but one cannot have both. 

6.4. Power spectrum results 

We apply the procedures described in Section 16.11 
to estimate the six possible CMB power spectra 
(TT,TE,EE,BB,TB,EB) from the QUaD maps. In ad- 
dition to the 100 and 150 GHz auto-spectra, we also estimate 
the 100-150 GHz cross-spectra as described in Paper II. In 
the case where the noise is uncorrelated between the two 
frequency channels, the cross spectra do not require the 
noise-debiasing step. Although, in practice, we do apply 
this correction, the correction is modest for the cross spectra 
so these measurements will be much less sensitive to the 
details of our noise model. The power spectrum results 
from PipeUne 1 are presented in Figure [SJ We make strong 
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Fig. 7. — Expectation values for the ACDM ii-mode power spectrum for 
both pipelines at 150 GHz. The mapping of true sky power to observed band 
powers for both pipelines is fully encoded in their associated band power 
window functions (Figure |6). 

detections of the TT, TE and EE spectra which show good 
agreement with those predicted by the best-fitting ACDM 
model to the 5-year WMAP results. The BB, TB and EB 
spectra are consistent with null within the noise. The results 
from Pipeline 2 show a similar agreement with the ACDM 
model. 

6.5. Jackknife tests 

We have subjected our analysis to the same set of systematic 
tests as performed in Paper 11. These tests involve splitting the 
data into two roughly equal parts. Maps made from the two 
data subsets are subtracted and the power spectra of the resid- 
ual maps are calculated. Any deviation of these "jackknife" 
spectra from null would indicate systematic contamination in 
the data. For a detailed description of the various ways in 
which we split the data, we refer the reader to Section 7 of 
Paper 11. 

The strongest test is the so-called "deck jackknife" test 
where we split the data according to the boresight rotation 
angle (deck angle) of the observations. This test, in particular, 
will be strongly sensitive to any residual ground contamina- 
tion remaining in the data after applying the procedure de- 
scribed in Section [3] In Figure |8] the deck jackknife spectra 
are plotted alongside the spectra measured from the undiffer- 
enced maps. It is clear from this figure that the power mea- 
sured in the deck-differenced maps is small compared to the 
measured signals in TT, TE and/oi EE. The other data splits 
which we consider are splitting the data according to orienta- 
tion of the PSB pairs on the focal plane (see Section|2]i, a split 
between the forward and backward scans and a split in obser- 
vation dates which roughly separates the data into the 2006 
and 2007 observing seasons. We also take the difference of 
the 100 and 150 GHz maps. This "frequency jackknife" test 
is not strictly a test for systematic effects in the data. Rather, it 
is a strong test for "foreground" (i.e. non-CMB) astrophysical 
emission. 

The cancellation of the signal apparent in Figure [8] is visu- 
ally impressive. To investigate whether the differenced spec- 
tra are formally consistent with null, as in Paper 11, we have 
performed t^sts against the null model. Note that we com- 
pare this statistic with the distribution as measured from 
our simulations rather than against the theoretical x^ curve. 



TABLE 1 

Jackknife PTE values from tests 



Jackknife 


TT 


TE 


EE 


BB 


TB 


EB 


Deck angle: 
100 GHz 
150 GHz 
Cross 


0.000 
0.008 
0.000 


0.415 
0.295 
0.028 


0.883 
0.963 
0.780 


0.933 
0.988 
0.197 


0.598 
0.258 
0.287 


0.917 
0.423 
0.527 


Scan direction: 
100 GHz 
150 GHz 
Cross 


0.008 
0.080 
0.000 


0.017 
0.665 
0.608 


0.122 
0.755 
0.155 


0.812 
0.153 
0.783 


0.478 
0.515 
0.487 


0.518 
0.485 
0.263 


Spht season: 
100 GHz 
150 GHz 
Cross 


0.743 
0.000 
0.273 


0.287 
0.387 
0.065 


0.350 
0.242 
0.110 


0.655 
0.022 
0.160 


0.840 
0.340 
0.630 


0.413 
0.647 
0.850 


Focal plane: 
100 GHz 
150 GHz 
Cross 


0.173 
0.530 
0.270 


0.872 
0.397 
0.012 


0.690 
0.910 
0.493 


0.813 
0.988 
0.105 


0.703 
0.933 
0.735 


0.672 
0.715 
0.578 


Frequency 
difference 


0.000 


0.362 


0.418 


0.588 


0.208 


0.783 



calculating the "probability to exceed" (PTE) the observed 
value by random chance. Low numbers therefore indicate a 
problem. The PTE values for each of our measured spectra 
(from Pipeline 1) are presented in Table[T] 

Examining the table, the PTE values for all the spectra bar 
TT reveal no significant problems, the numbers being con- 
sistent with a uniform distribution between zero and one. In 
contrast, many of our TT jackknife spectra are clearly incon- 
sistent with a null signal. The failure is perhaps excusable 
in the case of the frequency difference (since there are astro- 
physical reasons why this test might fail) but taken as a whole, 
the statistics for TT (and to a lesser extent, some of the TE 
numbers) suggest that there is some degree of residual sys- 
tematics present in our temperature maps. The PTE statistics 
for Pipeline 2, although not identical, show the same general 
pattern of jackknife failures for the TT spectra. Comparing 
to our previous results, there were hints of a problem with 
the TT jackknife tests in Paper 11 but to a lesser extent than 
is apparent now. This is possibly due to the fact that with 
our increased sensitivity we can measure both the signal and 
systematics to greater precision. Alternatively, it could indi- 
cate that the template removal procedure is not as effective as 
field-differencing in removing the ground. 

These jackknife failures indicate that residual systematic ef- 
fects in our temperature maps are significant with respect to 
the errors on the jackknife spectra. However, the jackknife er- 
rors contain only noise whereas the errors on our measured 
TT and TE CMB power spectra also contain considerable 
sample variance. To assess how significant the residual sys- 
tematics are, in Figure |9l we plot the measured TT and TE 
jackknife band powers alongside the signal band powers for 
the 150 GHz channel. Clearly, when compared to the sam- 
ple variance dominated signal spectra, the degree of residuals 
are much less significant. Repeating the x^ analysis using the 
signal CO variance matrices in place of the jackknife covari- 
ance matrices, we find that the residual contamination is neg- 
ligible. In summary, although our TT jackknife tests indicate 
the presence of residual systematics, they also clearly demon- 
strate that these residuals are irrelevant compared to both the 



Improved results from the QUaD CMB experiment 



11 




Fig. 8. — Full set of power spectrum results from Pipeline 1 w here the quantity pl otted in the y-direction is £(£+ l)Q/27r. The red curves are the spectra 
predicted by the best-fit ACDM model to the WMAP 5-year data (Dunkley et al."2009). The blue points s how the power spectra measured from our deck angle 
differenced maps and represent a stringent test for residual ground contamination — see the text of Section [631 for details. 
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TABLE 2 

X" TESTS OF THE COMBINED SPECTRA AGAINST 
THE ACDM AND NULL MODELS. 



1 ■ 





vs. ACDM 


vs. null 


f -^T J „ f . 

^ lor Lj a.o.i. . 






TT 
1 1 


1 S 71 

1 J. / 1 


1 SI 9 


1 c 




i iU. Jo 


EE 


24.87 


413.81 


BB 


37.00 


37.87 


TB 


28.99 




EB 


22.25 




PTE values: 






TT 


0.830 


< machine precision 


TE 


0.208 


7.9 X 10-'5 


EE 


0.303 


< machine precision 


BB 


0.024 


0.019 


TB 


0.145 




EB 


0.445 





measured sky signal and its associated sample variance. Note 
this is also true for the levels of foreground contamination in 
our maps since our frequency difference test is sensitive to 
foregrounds. 



Signal 
Deck angle 
Scan dirn. 
Split season 
Focal plane 
Frequency 



500 1000 1500 aooo 

multipole, I 

Fig. 9. — Signal and jackknife TT (top panel) and TE (bottom panel) power 
spectra for the 150 GHz channel. The black points show the measured signal 
spectra and the colored points are the jackknife spectra. The levels of power 
measured in the jackknife maps are negligible when compared to the sample 
variance driven errors on the signal spectra. 



6.6. Examination affinal power spectra 

To produce a set of final power spectra, we combine the 
single-frequency and 100-150 GHz cross spectra shown in 
Figure [8] following the procedure described in Section 9 of 
Paper 11. The final estimates of the power spectra are shown 
in Figure [To] for both of our analysis pipelines. Once again, 
the model curves plotted are the best-fitting ACDM model to 
the WMAP 5-year results. We have performed tests of all 
spectra and both analyses show perfectly acceptable agree- 
ment with this model and with each other Pipeline 1 and 
PTE values are given in Table |2l 

Comparing the two sets of results presented in FigurefTOlwe 
see that the nominal error-bars are smaller for Pipeline 2 and 



that the Pipeline 2 points appear to trace a slightly smoothed 
version of the Pipeline 1 points. Both of these effects are sim- 
ply a result of the di fferin g band power window functions as 
discussed in Section 16.31 Note that neighboring band pow- 
ers in Pipeline 1 are ~ 10% anti-correlated whereas they are 
~ 10% positively correlated in Pipeline 2. Correlations be- 
tween non-adjacent band powers are negligible for both anal- 
yses. 

Figure [TT] shows a comparison with our previous results 
from Paper 11. Note that we perform this comparison using 
Pipeline 2 since this pipeline is an extension of the analysis 
presented in Paper 11 and so is directly comparable. Two ef- 
fects are apparent in this figure. First, the uncertainties on 
all of our power spectra have been reduced by ^ 30% as a 
result of the increase in sky area afforded by our template- 
based ground removal technique. Second, implementing the 
improved beam models described in Section|4]has resulted in 
a slight increase in the amplitude of our power spectrum mea- 
surements for multipoles, £ > 700. This impacts mostly on 
the high signal-to-noise measurements of the TT spectrum on 
small scales. 

Figure[T2]shows our measured power spectra from Pipeline 
1 in comparison with the published results from a number of 
other CMB experiments. 

The QUaD power spectra data, along with the associated 
band power covariance matrices and band power window 
functions (for both of our analysis pipelines) are available for 
download at http://quad.uchicago.edu/quad 

6.7. Acoustic oscillations in the E-mode power spectrum 

We have assessed the significance with which we detect 
the acoustic oscillations in the EE power spectrum by re- 
peating the analysis of Section 9.3 in Paper 11. For this test, 
we compare our EE measurements against both the ACDM 
model and against a heavily smoothed version of the ACDM 
curve. The results of this test are shown in Figure [T3] The 
QUaD detection of acoustic oscillations in the E-mode power 
spectrum is now beyond question - the probability that the 
true E-mode spectrum is a smooth curve has dropped from 
0.001 with our previous analysis to < 10"'"*. We have also 
repeated our Paper 11 analysis where we used "toy mod- 
els" of the E-mode spectrum to fit the peak spacing, phase 
and amplitude of the acoustic oscillations. With our previ- 
ous measurements, we constrained the peak spacing to be 
A4 = 306 ± 10, the phase to be = 13° ± 33° and the am- 
plitude to be fl = 0.86 ±0.17. Repeating the analysis with 
our new measurements, we find A^j = 308 ± 7, = 6° ± 22° 
and a = 0.96 ± 0.10. For comparison, when we perform the 
analysis with the QUaD band power values replaced by their 
expectation values in the ACDM model, we find A£s = 310, 
(p = 13° and a = 0.99. These results confirm that the polar- 
ization peak spacing and the phase relationship between the 
temperature and polarization acoustic oscillations are as ex- 
pected in the ACDM model. Passing this (non-trivial) test 
further strengthens the foundations of this model. 

7. COSMOLOGICAL INTERPRETATION 

The QUaD results presented in this paper are the most pre- 
cise determination of the CMB polarization, and of its corre- 
lation with the CMB temperature, at angular scales £ > 200, to 
date. Within the standard ACDM model, given a precise mea- 
surement of the TT power spectrum (e.g. from WMAP), the 
TE and EE spectra on all but the very largest scales are de- 
terministically predicted. Nevertheless, sufficiently accurate 
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Fig. 10. — Final power spectrum results obtained from an optimal combination of the 100 GHz, 150 GHz and cross-frequency power spectra. Spectra from 
both pipehnes (again, plotted as l)Q/27r) are shown in comparison to the expected spectra in the concordance ACDM model. Note the TT power spectrum 
is plotted with a log-scale in the y-axis. For clarity, the two sets of points have been slightly offset in the horizontal direction. 



measurements of these spectra can still tighten constraints on 
cosmological parameters. In particular, precise measurements 
of the TE and EE acoustic peaks and valleys can help con- 
strain the cosmological parameters which determine this peak 
pattern. 

Beyond the standard ACDM model, power spectrum mea- 
surements from small-scale experiments such as QUaD can 
add further information. Two such extensions to the ACDM 
model which are well-motivated in the context of single-field 
slow-roll inflation models are a possible tensor component in 
the primordial perturbation fields and a scale dependence (or 
"running") in the scalar spectral index. Placing constraints 
on tensor modes through measurements of the large-scale 
Z?-mode polarization induced by a gravitational wave back- 
ground from inflation is a major goal of ongoing and future 
CMB polarization experiments. Although the scales probed 
by QUaD's polarization measurements are too small to con- 
strain this Z?-mode signal directly, QUaD can help through 
its ability to constrain the scalar spectra index, «„ since this 
parameter is correlated with the tensor-to-scalar ratio. For a 
ACDM model extended to include a running in ris, measure- 
ments of the high-£ temperature power spectrum can be useful 
in constraining both n, and the degree of running. 

In this section, we constrain cosmological models by 
adding the QUaD temperature and polarization data (i.e. our 
new measurements of the TT, EE, TE and BB spectra) to 
the results of two other CMB e xperiments — the WMAP 5- 
year analysis jNolta et al.ll2009l) and the final results from the 



ACBAR experiment dReichardt et al.ll2009l) . We will also in- 
vestigate the effect of adding large-scale structure data by in- 
cluding measurements of the present-day matter power spec- 
ti-um, P{k) fr om the SDSS L uminous Red Galaxy (LRG) 4th 
data release (iTegmark et al.ll200 6). 

In this paper, we focus on what the QUaD dataset taken 
as a whole adds to parameter constraints. In addition to 
investigating further extensions to the ACDM model, we 
will fully explore the consistency of our temperature-only 
and polarization-only parameter constraints in a future pa- 
per (Gupta et al., in prep). See also Castro et al. (2009) for 
temperature-only and polarization-only constraints obtained 
using our previous power spectrum results of Paper II. 

7.1. Methodology 

To obtain our constraints, we perform a Monte-Carlo 
Markov Chain (MCMC) sampling of the cosmological pa- 
rameter space. To do this, we use the publicly available Cos- 
moMC package (Lewis & Bridle 2002), which in turn uses 
the CAMB code (Lewis et al.ll2000() to generate the CMB and 
matter power spectra.^ 

^ Note that we have used the June 2008 version of the CAMB software. 
This version included a revised model of the reionization history as com- 
pared to previous versions of CAMB. In particular, the mapping between 
the optical depth, r and the reionization redshift, Zk changed at the ~ 10% 
level — see ILewij ^2008) for details. This should be borne in mind when 
comparing our results to previous analyses such as Dunklev et al. ( 2009) and 
iReichardt et al., t20Q^) who used pre-March 2008 versions of the CAMB soft- 
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We make use of the publicly available WMAP likelihood 
software from the LAMBDA** website. We marginalize over a 
possible contribution to the temperature power spectrum from 
the Sunyaev-Zel'dovich (SZ) effect following the WMAP 
analysis (Dunklev et al. 2009). To do this, we use the SZ tem- 
plates from Komatsu & Seliak(.200 2) (also available from the 
LAMBDA website) and the known frequency dependence of 
the SZ effect. In order to avoid possible contamination from 
residual point sources, we exclude the ACBAR band powers 
above £ = 2000. For the same reason, we do not include our 
own TT measurements at ^ > 2000 (Friedman et al. 206§). 

Marginalization over the WMAP beam uncertainty is in- 
cluded in the WMAP likelihood code and we also marginalize 
over the quoted ACBAR calibration and beam uncertainties. 
We take the latter to be a 2% error on a 5 arcmin (FWHM) 
Gaussian beam as assumed in the ACBAR CosmoMC data 
file. For QUaD, we marginalize over our 3.5% calibration un- 
certainty and over the uncertainty in our beam. As described 
in Appendix IaI our beam uncertainties are dominated by un- 
certainties in the level of our sidelobes rather than in the effec- 
tive FWHM of our main lobe beams. We therefore marginal- 
ize over the full ^-dependent beam uncertainty shown in Fig- 
ure[T9] Where we include the SDSS LRG data, we marginal- 
ize over both the amplitude of the matter power spectrum and 
over a correction for scale-dependent non-linear density evo- 

ware. 

* http://lambda.gsfc.nasa.gov/ 



lution using the methods described in lTegmark et a l.' ('20061). 

We model the likelihood functions fo r the QU aD auto- 
spectra as offset log-normal distributions ("B ond et al.ii2000l) . 
The required noise-offsets are derived from our signal-only 
and noise-only simulations. (We model the TE likelihood as 
a Gaussian distribution.) We include all covariances appar- 
ent (above the numerical noise) in our simulation-derived co- 
variance matrix (equation [T6]). In addition to same-spectrum 
covariances, this includes non-zero TT-TE, EE-TE, TT-EE 
and EE-BB correlations. 

Note that, for our main MCMC analysis, we do not include 
our measurements of the parity-violating spectra, TB and EB 
since these spectra are expected to vanish in standard ACDM 
models and its usual variants. However, in Section 17.71 we 
will use these spectra to constrain possible parity-violat ing in- 
teractions to the surface of last scatt ering (se e e.g. L ue et al.l 
[T999I) following our previous work dWu et al.l20()9l) . 

Finally, in Section 17.81 we use our polarization measure- 
ments to place a formal upper limit on the strength of the 
lensing B-mode signal. 

Our basic ACDM cosmological model is char- 
acterized by the following six parameters (where 
/i = //o/ [lOOkms"'Mpc"'], Hq being Hubble's constant 
in units, kms~'Mpc~'): the physical baryon density, 57;,/?^; 
the physical cold dark matter density, VLji^; the ratio of 
the sound horizon to the angular diameter distance at last 
scattering, 9 = VsIDa, the optical depth to last scattering. 
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Fig. 12. — QUaD measurements of the TT, TE EE and BB power spectra compared to results from the WMAP (Nolta et al."2009), ACBAR fReichardt et alj 
|2009), BICEP (Chiang et al. 2009), B03 iPiacentmi et al. 2006; Montrov et al 2006), CBI iSievers et al. 2007), CAPMAP ( Bischoff et al. 2008), MAXIPOL 
TWu et al..2007.) and DASI iLeitch et aL2005.) experiments. The BB measurements are plotted as 95% upper limits. The smooth black curves in each panel are 
the power spectra expected in the best-fit ACDM model to the WMAP 5-year data. 

and which are derived from this basic set, are the dark energy 
density, (assumed here to be a simple cosmological 
constant), the age of the universe, the total matter density, 
Vl,„, the amplitude of matter fluctuations in 8 /i~'Mpc spheres, 
(78, the redshift to reionization, Zre and the value of the 
present day Hubble constant, Hq. For all our analyses, 
we assume a flat universe and include the effects of weak 
gravitational lensing. We impose the following broad priors 
on our base MCMC pai-ameters: 0.005 < ^Ibh' < 0.100; 
0.01 < n^h^ < 0.99; 0.5 < 6 < 10.0; 2.7 < A < 4.0; 
0.5 < «s < 1.5; 0.01 < T < 0.80. There is also a prior 
imposed on the age of the universe (10 < Age [Gyrs] < 20) 

and on the Hubble constant (40 < Ho [kms"'Mpc"'] < 100). 

We also investigate models extended to include both a run- 
ning in the scalar spectral index, = dtis/dlnk and/or a 
possible tensor contribution. Assuming a power law for the 
tensor modes, V, oc k"' , we parametrize their amplitude by the 
tensor-to-scalar ratio, r = Vt/Vs- We adopt a uniform prior 
measure for r between and 1 . For the running spectral index 
model, we adopt a prior of -0.5 < Wrun < 0.5 on the running. 
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Fig. 13. — QUaD's measurements of the EE spectrum (black points) com- 
pared to the ACDM model (red curve) and a model without peaks (green 
curve). The data are incompatible with the no-peak scenario — the probabil- 
ity that the smooth curve is connect is < 10"'''. 

r; the scalar spectral index, Wj; and the scalar amplitude. 
As = In [lO'°Aj] . Here, A, is the amplitude of the power 
spectrum of primordial scalar perturbations, parametrized by 
Vsik) =As{k/kl)"~^. We discuss the choice of pivot-point, kl, 
in the following section. Other parameters which we quote. 



7.2. Choice of scales ( "pivot-points" ) for presentation of 

results 

For the primordial power spectrum parametrization which 
we have chosen, we need also to choose a scalar pivot-point. 
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TABLE 3 
Basic 6 Parameter Constraints 



WMAP 



WMAP+ACBAR WMAP+QUaD WMAP+ACBAR+QUaD WMAP+ACBAR+QUaD+SDSS 



A 

Qa 
Age 

Zre 



n rms+0 0006 
1.0406^'®^o 



().090^„„ 
965+<' '"'* 

72.1+2« 



Q 111+0.006 

1 . n nn27 



0090^.017 
0.966^:«'3 
3 12+0.04 
74+8:8i 

"■'^-0.03 

13 63+0.12 

0.26-; 
0.80|i 

10.5+['^ 

+2:5 
-2.5 



0.0227+ 



1.0403 



72.3: 



7+0.0005 
'-0.0005 

0-090^.oif 
962+<""3 

3i;^.-8i3 

13 69+0.12 

25+o"3 
-Q.!!3 

"•'^-0.03 
10.5+1-* 



72.5: 



0227+<'-'""'5 
Q lQg+O.005 

1^0415^52 



0.089^0,7 
962+<' '"3 

3 11+0.04 

75+^:8! 
'-'•'■'-0.03 
1 3 66+* ' 1 

25^-'^^ 
-Q.Q3 

10.5+'-^ 



72.4- 



n rm7+0.0005 
O.lOS+flOOl 



1 n414+<'"''22 

962+" '"2 
3,lM^ 

-fl.Q3 
75+0 02 
'-'-0.02 

13 66+0- 10 
10.00 

25+°"2 
-Q.Q2 
79+0.02 
'=^-0.02 

72.7+{;^ 



Note. — We quote the scalar amplitude as A., = In [lO'^A.,] for a pivot-point of kl = 0.013 Mpc 



kl, the wavenumber at which Wj and A, are evaluated. Within 
standard ACDM, Hs is modeled as independent of scale and 
we can map constraints on A, obtained at one pivot-point to 
an arbitrary new pivot-point, kl', using 

A,(ki')=AAIcl)il<'/0"'~'- (21) 

For models including a running in the spectral index, both 
n, and A, are dependent on scale. For these models, we can 
map constraints from an old to a new pivot-point using 

AAK')=AAK) (;t^7;t^)«««)-l+l"™l"«'A:)^ (22) 

nAK') = nAK) + n,u„ 111(^:7^:) • (23) 

Correlations between the two parameters, and Wnin 
are dependent on the pivot-point at which one chooses to 
present results. In particular, there is a scale at which 
the un certainties on these two parameters become uncorre- 
cted (Copeland et al.l [19981) . Choosing to present results at 
this "decorrelation scale" has the attractive feature that the 
marginalized ID constraint on i s not degrad ed by allow- 
ing the running to be non-zero. iFinelli et al.l (12006 ) pre- 
sented parameter constraints from CMB and large-scale struc- 
ture da ta using a pivo t-point of k^ = 0.01 Mpc~' whereas 
iPeiris & Easther ( l2006h identified a decorrelation scale of 
ki w 0.02 Mpc"' using the WMAP 3-yr data. 

In order to find th e decorrelatio n pivot-point, we have fol- 
lowed the analysis of lCortes et al.l (j2007) who describe a tech- 
nique to fit MCMC chains for the decorrelation scale. They 
found a decorrelation scale of kl = 0.017 Mpc"' using the 
WMAP 3-yr data set. Repeating their analysis using the 
WMAP 5-yr data, we find kl = 0.013 Mpc"'. For simplic- 
ity, we choose to present our constraints on As and ris at this 
decorrelation scale for all of the models and dataset combina- 
tions which we have investigated. 

For models including a possible tensor component, we still 
quote our constraints on A,, and at kl = 0.013 Mpc"' but 
for the tensor-to-scalar ratio, r, we use a tensor pivot-point of 
k!^ = 0.002 Mpc"' . We do not attempt to remap our constraints 
on r to a more optimal pivot-point since the only meaningful 
data contributing to a constraint on r is the WMAP tempera- 
ture power spectrum on very large scales (for which, a tensor 
pivot-point of kl = 0.002 Mpc"' is appropriate). Note also that 
for these models, we enforce both the first and second infla- 
tion consistency equations (e.g. ,Lidsev et al...l997.) : r = -8n, 



and = n, [n, - (wj- 1)]. Additionally enforcing the second 
equation ensures that the first con sistency equation holds to 
linear order in A In A: on all scales dCortes et al.ll2007h . 

7.3. The concordance ACDM model 

In Table [3] for each dataset combination, we list the mean 
recovered values for each parameter, along with their associ- 
ated 68% confidence limits, marginalized over all other pa- 
rameters. In Figure [14] we plot the marginalized ID con- 
straints for the WMAP, WMAP + QUaD and WMAP + 
ACBAR combinations. Cleai'ly, the WMAP data dominates 
when we add in either the ACB AR or the QUaD d ata, as was 
found in our previous analysis dCastro et al.ll2009l) . However, 
the addition of either of these data sets does provide additional 
information on fl^h^, Vlch^ and 6. The biggest improvement in 
constraints is in where the WMAP + ACBAR + QUaD com- 
bination tightens the limits by ^25% compared to WMAP 
alone. This additional constraining power comes mostly from 
the QUaD data. 

The mean values of these parameters also shift a little but 
the only significant discrepancy is perhaps in the recovered 
value of 6. Here, we find the WMAP + ACBAR combi- 
nation prefers a somewhat hig her value (in ag reement with 
ACBAR's own analysis; Reic hardt et al.l [20091) whereas the 
addition of QUaD data does not change the WMAP-only pre- 
ferred mean value but simply tightens the constraint. 

Note that in comparison to the WMAP team's own analy- 
sis (Dunkley et al. 2009), we recover slightly different mean 
values for r and more significantly different values for Zre- 
This is due to the different reionization model used in the 
later version of the CAMB software which we have used. We 
note in passing that the majority of the constraining power 
in the QUaD data comes from the measurements of the po- 
larization po wer spectra as found with our previous analysis 
dCastro etal.li2009) . 

7.4. Running spectral index model 

The ID and 2D marginalized constraints on our base pa- 
rameters for the running spectral index model, as obtained 
from our WMAP + QUaD runs are shown in Figure [15] along 
with the constraints using only the WMAP data. The recov- 
ered parameter values and their uncertainties are listed in Ta- 
ble[4] 

The impact of QUaD data is greater for this model — 
the QUaD data adds significantly to the constraints on fii/i^. 
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Fig. 14. — ID likelihood distribution for our base MCMC parameters for 
the basic 6-parameter ACDM model. The constraint on the scalar amplitude 
is presented at a pivot-point of kl = 0.013 Mpc"' . 

rieh^, 6, and Wmn, reducing the marginalized ID errors on 
these parameters by up to 20%. Adding both the QUaD and 
ACBAR data has an even greater impact, reducing the errors 
on these parameters by up to a third compared to the WMAP- 
only uncertainties. Of particular interest are the constraints 
in the n, - «iun plane since many theories of inflation predict 
both a deviation from n, = 1 and/or a small negative running. 
Constraints from the WMAP + ACBAR + QUaD combina- 
tion are shown in the left hand panel of Figure [16] together 
with the constraints from WMAP alone. Our ID marginal- 
ized constraint on the running from the combined data set 
is Mrun = 0.046 ±0.021, 2.2(7 away from the Wi-un = model. 
Adding the LSS data to the mix improves the constraints even 
further, in our analysis tightening the Icr error on n^un from 
0.021 to 0.018. The significance of a non-zero running is also 
reduced on addition of the LSS data. 

Comparing Tables |3] and H] we see also that the ID 
marginalized constraint on the spectral index, «j is not weak- 
ened by allowing a non-zero running. This is due to our u se of 
the decorrelation pivot-scale as described in Section IT2l The 
results also show that the constraints on Hs obtained for the 
standard 6-parameter ACDM model are robust to marginal- 
ization over a possible running. For example, with the WMAP 
+ ACBAR + QUaD combination, the constraint on goes 
from «, = 0.962 ± 0.013 to «, = 0.965 ± 0.013 when we al- 
low for a possible non-zero running. For comparison, if in- 
stead, we use the WMAP-preferred pivot-point (kl = 0.002 
Mpc~'), the marginalized ID constraint on n, is degraded to 
«s = 0.962 ± 0.019 in the presence of running. 

Hints of a negative running in the spectral index have been 
observed in previous CMB analyses (e.g. Dunkley et al. 2009; 
iReichardt et al. 2009). With the addition of the new QUaD 
data, this suggestion of a negative running not only persists, 
but is strengthened. We do, however, stress that the combined 
result shown in the left hand panel of Fig.[T6]is still consistent 
with zero running at the 3a level. Nevertheless, it is worth 



examining the implications for inflation models if the running 
was as large and as negative as the best-fit value returned from 
our analysis of the combined CMB data set. In this respect, 
iMalquarti et al. (2004) have pointed out that an observed run- 
ning of «iun < -0.02 would effectively rule out large field in- 
flation models. More generally. lEasther & PeirisI (12006 ) have 
demonstrated that for a large negative running, single field 
slow roll inflation models will last less than 30 e-folds af- 
ter entering the horizon. This amount of inflation is insuffi- 
cient if inflation happened at the GUT scale. Consequently, 
an observation of Wrun ~ -0.05 would require inflation the- 
ory to move beyond the simplest models, e.g. by considering 
multiple fields and/or modifications to the slow-roll formal- 
ism (e .g. IChung et alj|2003l: lMakarovll200l iBallesteros etaLl 
l2006h . 

7.5. Tensor modes 

Our constraints for the ACDM model including a possible 
tensor component are listed in Table |5] in terms of the mean 
recovered parameter values and their uncertainties. For r, we 
quote the 95% one-tail upper limit since as expected, no de- 
tection of tensors is made. For our WMAP-only analysis, we 
recover a slightly weaker limit (r < 0.4 8) than that obtained b y 
the WMAP team themselves (r < 0.43: lDunklev et al.l2009h .^ 
Adding either the ACBAR or QUaD data, this is reduced to 
r < 0.40. 

The WMAP + ACBAR + QUaD combination produces a 
constraint on tensor modes of r < 0.33, the strongest from 
the CMB alone to date. Note that this constraint does not 
come from our upper limits on the BB spectrum. It is, in fact, 
driven by a preference of the small-scale data (particularly the 
QUaD EE and TE data) for a somewhat lower spectral index 
compared to that preferred by WMAP alone — a lower «, al- 
lows more of the large-scale TT power observed by WMAP 
to come from scalar perturbations and therefore the maximum 
allowed tensor contribution is reduced. Our CMB-only con- 
straints in the r-W j plane are plotted in the right-hand panel of 
Figure [16] 

7.6. Running spectral index and tensor modes 

When we allow for both a running in the spectral index and 
a tensor contribution the constraints weaken considerably ver- 
sus either on their own. In the left-hand panel of Figure [17] 
for the WMAP + ACBAR + QUaD combination, we plot con- 
straints in the «iun-n.! plane with and without marginalization 
over a possible tensor component. The right panel of this 
figure shows the corresponding constraints in the r-n^ plane 
with and without marginalization over a possible running in 
the spectral index. For this model, the addition of QUaD 
and ACBAR data still improves the constraints in the spec- 
tral index running (and indeed, still strongly suggests a small 
negative running) but the constraints in the r-n^ plane in the 
presence of running do not improve on the WMAP-only re- 
sult. This degradation in the constraints on r when we allow 
for a running in the spectral index is further demonstrated in 
Table [6] where we quote the ID marginalized constraints on 
the parameters {r^n^, nmn} for the tensors-only, running-only 
and tensors + running models. In this table, we present the 
results for the WMAP + ACBAR + QUaD combination and 
for the case where we add in the SDSS LRG data. 

' Repeating our MCMC analysis using the pre-March 2008 version of 
CAMB and adopting WMAP's choice of both scalar and tensor pivot-points, 
we recover a result consistent with the WMAP analysis. 
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TABLE 4 

Parameter Constraints for the running spectral index model 



WMAP WMAP+ACBAR WMAP+QUaD WMAP+ACBAR+QUaD WMAP+ACBAR+QUaD+SDSS 

n,h- 0221+* ™'' 0221+0 0™ 0219+0 ™' n 0219+<' ''<""' 0223+" '""'* 



n,/j2 o.iie^o™ 0.120^™ 0.117^™ o.i20!<!™ o.iii 



Hq 68.9!|| 68.0!|| 68.5!^| 67.5!|^ 



^ ^ ^ ^ _ ^ +0.004 

'^0 1 (14 14+^'^^^' 1 m99+^'-^^^'* 1 n409+^-^^^3 1 04n+^''^'^^2 
^-"^^ -0.0027 '•"-'"-0.0024 '-"^"^-0.0023 '■"^'-'-0.0022 

0.964jf, 0.9671^1 0.963^«jJ 0.965^«j3 oM^^ 

-0.031^gji -0.040^1 -0.038!Sjj| -0-046!gJl -0.028!S]jl« 

J o 1 £+0.06 o 1 o+fl.05 T 1 7+0.06 T 1 n+O.05 a i <;+fl.04 

■'■•"-0.06 -'■•°-0.05 ■^■•'-0.06 ■^••^-0.05 -'■•■'-0.04 

Ha 0.70^^«^ o.69!«j^ 0.70^«^ q gg^^os o.74:"^02 

Age 13.81!!>J 13.78!ijJ 13.82^| 13.82!i>| 13.72^| 

n,„ o.3ot!!^ o.3i;U;J!i o.3ot!!§ 0.3 

1-005 "-26^002 

OM^ot 0.83;«;OJ 0.81^^^ 0.83t| 0.79+»:0| 

z„ 11.2+|'^ 11.4+1-f 11.4+!-^ 11.7+1-f 11.2+:' = 



Note. — The pivot point used for As and Wj is /cj = 0.013 Mpc 



7.7. Constraints on parity violation 

In the preceding sections, we used the QUaD TT.EE ,TE 
and BB spectra to constrain the parameters of standard ACDM 
models and its usual extensions. For that analysis, we did 
not use our measurements of the cross-polarization spectrum 
(EB) or the correlation of temperature with B-modes (TB), 
since these spectra are expected to vanish in a universe which 
respects parity conservation (which the above models do). In 
this section, we use these two spectra (along with the TE, EE 
and BB measurements) to constrain a possible parity violation 
signal on cosmological scales. In the presence of parity vio- 
lating interactions, a rotation in the polarization direction of 
CMB photons will be induced as they propagate from the sur- 
face of last scattering. If parity violating effects are present 
on cosmological scales, there will therefore be a net local 
rotation of the observed Stokes parameters, Q and U in the 
measured polarization map. This will mix E and B modes 
resulting in non-zero expectations for the TB and EB power 
spectra. Parametrizing the parity violation effect with a ro- 
tation angle, Aa, the expectation values for the TB and EB 
spectra in terms of the cosmological TE and EE spectra are 
given by: 

C[* = C[^sin(2Aa) (24) 
Cf^ = lcf^sin(4Aa). (25) 

In addition, assuming that primordial and lensed B-modes are 
negligible (which is an excellent assumption given our sensi- 
tivity), the expectation value for the BB spectrum is 

C^'^ = Cpsm^(2Aa). (26) 

We used our previous results to place a constraint of Aa = 
0.53° ±0.82° ±0.50° dWu et aljliool where the two quoted 
uncertainties are the random and systematic components re- 
spectively. 

Here, we repeat this analysis with our new measurements. 
The analysis is model-independent in the sense that we con- 
struct our estimator for Aa in terms of the observed power 

The pai'ity violation effect is completely degenerate with an error in 
the calibrati on of the polariza tion co-ordinate system of the experiment. As 
described in IWu et al.l {20091) . we are confident in our calibration to at least 
0.5°. 



spectra and do not assume a cosmological model for the EE or 
TE signals. For details of our estimator and analysis (which 
has jiot changed since our previous work), we refer the reader 
to iWuetal.l(i2()0l . 

We apply the estimator to the real QUaD data and assign 
error-bars due to random noise and sample variance by pro- 
cessing the suite of simulations containing both signal and 
noise through the analysis. The result is shown in Figure [TS] 
where we plot both the results from simulations and from the 
real data. Note that our simulations contain no parity violat- 
ing signals and so should scatter about zero, which they do. 
We take the scatter in the results from the simulations as our 
random error. Adding in the systematic error, our final result 
is 

Aa = 0.64 ± 0.50 (random) ± 0.50 (systematic) . (27) 

The random error has been reduced by ^ 40% with respect to 
our previous analysis in line with expectations. 

Our result can be compared to the li mits obtained from the 
WMAP 5-year data (Aa = -1 .7 ± 2. 1 : iKomatsu et al.ll2009h 
or to the limits obtained from the combination of th e WMAP 
5-yea r data and the B03 results (Aa = -2.6 ± 1.9; IXia et all 
120081) . We note that both of these quoted results include ran- 
dom errors only and do not include estimates of the systematic 
errors on the WMAP and B03 polarization calibration angles. 
Even when we include this systematic uncertainty for QUaD, 
our result is clearly a marked improvement over these previ- 
ous analyses. 

7.8. Limits on the lensed B-mode signal 

Although QUaD has not made any detection of B-modes, 
it is the most sensitive small-scale CMB polarization experi- 
ment to date. We can therefore place the leading upper limit 
on the presence of a small scale B-mode signal. The signal 
expected to dominate on the scales at which QUaD is sensi- 
tive is that induced by gravitational lensing of E-modes by 
intervening large scale structure. As well as measuring cos- 
mological B-modes from inflation, future polarization experi- 
ments will target this lensing signal from which useful infor- 
mati on can be gained on dark energy and ma ssive neutrinos 
(e.g. lKaplinghatet al. 2003; S mith et al.ll2006l) . 

Assuming a single flat band power between ^ = 200 and 
I = 2000, we find £(i+ \)Cf /It: = 0.17 ± 0.17 fjK- with a 
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Fig. 15. — ID and 2D likelihood distributions recovered for the base MCMC parameters for the ACDM model extended to include a possible running in the 
scalar spectral index. The constraints on n, and A, are shown for the decorrelation pivot-point of kl = 0.013 Mpc"' . In the 2D panels, we indicate the regions of 
parameter space which enclose 68% and 95% of the likelihood as the inner and outer contours respectively. The results for the WMAP + QUaD combination are 
shown over-plotted on the WMAP-only results. Adding the QUaD data tightens the constraints on Q,th^, flch^, «iun and 9 by up to 20%. 



95% upper limit of 0.57 /iK^. The errors quoted are estimated 
from the scatter in the results obtained from the suite of simu- 
lations containing both signal and noise. For comparison, the 
ACDM expectation value for this band power is 0.058 /iK^. 
Alternatively, assuming the ACDM shape for the lensing sig- 
nal, and simply fitting for its amplitude between £ = 200 and 
£ = 2000, our constraint on the amplitude" is 2.5 ±4.5 with a 
95% upper limit of 12.5. 

Note that although we have used all of our BB band powers 
to obtain the above constraints, the window function of our 
estimator is strongly skewed towards lower multipoles where 
the band power uncertainties are much smaller The effective 
range of multipoles to which our upper limits apply is, in fact, 

' ' Our normalization convention is such that the amplitude of the lensed 
B-mode signal in the concordance ACDM model is unity. 



170 < £ < 400 rather than the nominal 200 < ^ < 2000 range. 

Although our 2a upper limits are an order of magnitude 
larger than the expected ACDM signal, they are, in turn, 
roughly an order of magnitude better than previously reported 
limits on the amplitude of the Z?-mode signal in this angular 
scale range. 

8. CONCLUSIONS 

We have presented a re-analysis of the final dataset from the 
QUaD experiment, a CMB polarimeter which observed the 
CMB at 100 and 150 GHz from the South Pole between 2005 
and 2007. A major part of this re-analysis was the develop- 
ment of a new technique for removing ground contamination 
from the data. The ground signal seen in QUaD data is polar- 
ized and, if not removed, contaminates all of the CMB power 
spectrum measurements. Our new procedure, which is based 



20 



QUaD collaboration - M. L. Brown et al. 



TABLE 5 

Parameter constraints including a possible tensor component 
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Note. — The pivot point used for As is kl 
ytl = 0.002 Mpc"'. 



: 0.013 Mpc while the pivot point used for the tensor-to-scalar ratio, r is 





t With tensors 
t Without tensors 




t With running 
I Without running 




Fig. 16. — Left panel: 68% and 95% confidence regions in the n.5-«run 
plane, marginalised over all other parameters, for the WMAP + ACBAR -I- 
QUaD combination as compared to those obtained from WMAP alone. The 
constraints are shown for a pivot-point of kl = 0.013 Mpc"' . No tensor com- 
ponent was allowed for either set of constraints. The constraints tighten by 
about one third. The mean recovered values also shift further away from the 
simple {«.5,«iim} = {1,0} model. Right panel: Marginalized constraints on 
the inflation parameters, r and n, from WMAP data alone and adding in the 
ACBAR and QUaD datasets. No running in the spectral index was allowed 
for these fits and the tensor-to-scalar ratio, r is presented for a tensor pivot- 
point of k[ = 0.002 Mpc"' . Once again, the inner and outer contours indicate 
the regions of parameter space enclosing 68% and 95% of the likelhood re- 
spectively. The 95% upper limit on r is reduced from r < 0.48 to r < 0.33. 
This constraint is driven by the preference of the additional datasets for a 
lower spectral index than is recovered from the WMAP data on its own. 



TABLE 6 

Constraints on inflationary parameters 



Parameter 



Tensors 



Running Tensors + Running 



CMB only: 

r 

dris/dink 



< 0.33 (95% c.l.) 



0.978+: 



-0.046:*;02| 
0.965^013 



< 0.60 (95% c.l.) 
-0.063^-055 

0.997^-026 



CMB + LSS: 

r 

dris/dink 



< 0.27 (95% c.l.) 



0.973+ 



-0.028+^ 



0.967+ 



< 0.61 (95% c.l.) 
-0 052+O-023 



Fig. 17. — Left panel: Constraints are shown (as 68% and 95% confidence 
regions) in the «run-'i.i plane with and without marginalization over a possible 
tensor component. These fits are for the WMAP -I- ACBAR -I- QUaD com- 
bination. Allowing a non-zero tensor component weakens the constraints 
considerably. However the addition of both QUaD and ACBAR to WMAP 
still favors a small negative running. Right panel: The 68% and 95% central 
confidence regions in the r-Us plane (for WMAP + ACBAR -f QUaD) with 
and without inarginalization over a possible running in the spectral index. 
Allowing the spectral index to run degrades the constraints to such an extent 
that the addition of QUaD and/or ACBAR data yields no improvement over 
the WMAP-only constraints. «j is evaluated at kl = 0.013 Mpc"' (for both 
panels) and r is evaluated at kl = 0.002 Mpc"' . 

on constructing, and subsequently subtracting, templates of 
the ground signal has allowed us to reconstruct maps of the T, 
Q and U Stokes parameters over the full sky area. Although 
the method is not entirely lossless, it provides, on average, a 
30% increase in the precision of the power spectra compared 
to our previous analysis which used field-differencing to re- 
move the ground. 

Through further detailed analysis of calibration data, we 
have also significantly improved our understanding of the 
QUaD beams. We have implemented new beam models 
which explicitly incorporate the effects of sidelobes, resulting 
in an increase of ~ 10% in the amplitude of our power spec- 
tra measurements for multipoles, I > 700. The shift in power 
is most relevant for our high-^ temperature power spectrum 
measurements where the signal-to-noise is high. 

We have presented results using our two independent anal- 
ysis pipelines. Though there are significant differences in the 
approach between the two pipelines, the final results agree 
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Fig. 18. — Constraints on possible parity violation interactions on cosmo- 
logical scales, parametrized in terms of the parity violation rotation angle, 
Aa. The histogram shows the estimates of Aa as measured from our suite 
of signal + noise simulations. The vertical red line shows Act as measured 
from the real QUaD data. The full and dashed blue lines show the 68% and 
95% confidence regions about the central value as estimated from the scatter 
in the simulation results. 

very well. Testing the power spectra against the best-fit 
ACDM model to the WMAP 5-year data, we find good agree- 
ment. Our measurements of the E-mode polarization spec- 
trum, and of the cross-correlation between the E-modes and 
the CMB temperature field, are the most precise at multipoles 
£ > 200 to date. Our measurement of the temperature power 
spectrum at ^ > 1000 is among the best constraints on temper- 
ature anisotropics on small angular scales and is co mpetitive 
with the final ACBAR result (Reic hardt et al.ll2009h . 

We have subjected our results to the same set of rigorous 
jackknife tests for s ystematic effects as was performed in our 
previous analysis dPrvke et al.l 12009 ). We find no evidence 
for residual systematic effects in our polarization maps. Al- 
though formally, many of our TT jackknife tests fail, the in- 
feiTed levels of residual systematics are negligible compared 
to our sample-variance driven error bars. Moreover, the very 
small level of power seen in our frequency difference maps 
and power spectra indicate that foreground contamination is 
also negligible compared to our uncertainties. 

We have used our power spectra measurements, in combi- 
nation with the WMAP 5-year results and the ACBAR results 
to place constraints on the parameters of cosmological mod- 
els. For the standard 6-parameter ACDM model, the QUaD 
data adds only marginally to the constraints obtained from the 
WMAP data alone. The impact of the QUaD data is greater 
in a model extended to include a running in the spectral in- 
dex, reducing the uncertainties in ilhh^, flch^, and ni-un by 
up to 20%. The addition of both QUaD and ACBAR data is 
more powerful still, improving the constraints on these four 
parameters by up to one third. For a ACDM model extended 



to include a possible tensor component, we find that the addi- 
tion of both ACBAR and QUaD data reduces the upper limit 
on the tensor-to-scalar ratio from r < 0.48 to r < 0.33 (95% 
C.I.). This is the strongest limit to date on tensors from the 
CMB alone. The improvement is driven by a tendency of the 
QUaD data to prefer a somewhat smaller spectral index than 
is infeiTed from WMAP data alone. 

We have used our measurements of the TB and EB power 
spectra to put constraints on possible parity-violating interac- 
tions on cosmological scales. Following our previous analysis 
(|Wu et al. 2009), we constrain the rotation angle due to such a 
possible "cosmological birefringence" to be 0.64° ± 0.50° ± 
0.50° where the errors quoted are the random and system- 
atic contributions. Our result is equivalent to a constraint on 
isotropic Lorentz-violating interactions of A:|^*)Qg < 1.5 x 10"^-' 
GeV (68% C.I.). 

Finally, we have placed an upper limit on the strength of 
the lensing B-mode signal using our measurements of the BB 
power spectrum. Assuming the concordance ACDM shape 
for lensing Z?-modes, we constrain its amplitude (where the 
normalization is such that the ACDM model has amplitude 
= 1) to be 2.5 ±4.5 with a 95% upper Hmit of 12.5. Alter- 
natively, assuming a single flat band power for £ > 200 we 
find a 95% upper Hmit of £(i+ \)Cf /It: < 0.51 iiK? for the 
amplitude of B-modes. 
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APPENDIX 



A. BEAM UNCERTAINTIES 

As described in Section |4] our new beam models involve either fitting the QUaD physical optics (PO) beam models to QSO 
data (Pipeline 1) or measuring the sidelobes directly from QSO maps under the assumption that the sidelobes are azimuthally 
symmetric (Pipeline 2). Although the predicted radially averaged beam profiles from both of these approaches appear to match 
the data very well, the fits are not perfect and are subject to an uncertainty in the sidelobe levels. There is also an uncertainty on 
the width of the main lobe, dominated by small temperature-dependent variations. Based on the fluctuations in the beam widths 
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multipole, I 

Fig. 19. — Fractional uncertainties on our beam transfer functions. The dashed curves show the errors due to the 2.5% uncertainty in the width of the main lobe 
beams. The dotted curves show the sidelobe uncertainties and the full curves show the total uncertainties. The con'esponding curves for the combined spectra 
are similar to the 150 GHz curves shown in red. At 150 GHz, the uncertainty in the level of the sidelobes dominates the beam uncertainty for the full £-range 
presented in this paper 



seen in our "rowcal" data'^, we estimate the remaining uncertainty on the main lobe width to be 2.5% of the effective FWHMs of 
5.2 and 3.8 arcmin at 100 GHz and 150 GHz respectively. We obtain the uncertainty on the level of our sidelobes from the errors 

returned from fitting our PO simulations to the QSO observations in Pipeline 1 . 

To propagate these errors onto uncertainties in the transfer functions of Section l673l for the error in the main lobe, we simply 
note that the effect of a fractional error, 5 in the FWHM of a Gaussian beam is well approximated by 

= exp [al(S^ + 2S)e(£+ 1)] - 1 , (Al) 



where ai, = ^fwhm / \/8 In 2 is the beam width. For the errors in the sidelobes, we take the minimum and maximum sidelobe levels 
as returned from the fits of the PO models to the data, coadd the resulting beam models across detectors and radially average to 
produce the minimum and maximum allowed radial profiles, B(9), for each frequency. Taking the Legendre transform of these 
profiles, 

Be = 2n J B(e)Pecos(9)dcos(e) , (A2) 
we estimate the error in our beam transfer functions due to the uncertainty in the sidelobes as 

= ^muxj " -^Ln/ • (A3) 

We take the quadrature sum of the errors due to the main lobe and sidelobe uncertainties to be the final error. These uncertainties 
are shown in Figure [19] along with the quadrature sum. Since our combined spectra are dominated by the 150 GHz channel, the 
curves for the combined spectra are approximately the same as the 150 GHz curves. 

B. ABSOLUTE CALIBRATION UNCERTAINTY 

As described in Section |5] our calibration is performed by taking the ratio of cross-spectra between the QUaD and B03 
temperature maps. This process is subject to several uncertainties. 

Firstly, there is a statistical error in the calibration ratio, predominantly due to noise in the B03 maps. To estimate this error, 
we perform Monte-Carlo simulations of the absolute calibration process. Assuming white noise and a ACDM power spectrum, 
we use the B03 hit-maps along with their stated sensitivities to produce simulations of the B03 maps. We do the same thing for 
QUaD and apply a known fiK V calibration to the simulated QUaD maps. Each pair of simulated maps is then passed through 
the absolute calibration analysis. The scatter in the calibration factors recovered from these simulations is 0.6% and we take this 
as the statistical uncertainty in our calibration. 

Our calibration ratio as a function of multipole is not perfectly flat but fluctuates about a mean value. Although some of this 
scatter will be due to noise (which is included in our estimate of the statistical eiTor), we conservatively also include this scatter, 
which we measure to be 1.1%, in our e rror budget. In addition, we have also performed the calibration analysis using each of 
the jackknife splits described in Section l63] Although the scatter found in the recovered calibration numbers do not indicate any 
significant inconsistencies, we also include this scatter in our error budget. 

These calibration data consisted of scanning each row of pixels in the focal plane across the bright HII region, RCW38 and were taken daily throughout the 
QUaD observations. Although RCW38 is not a true point source, the fluctuations in the per-channel beam widths put a tight constraint on temperature dependent 
seasonal fluctuations in our main lobe beams. 
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TABLE 7 

Absolute calibration uncertainties for QUaD 



Source Uncertainty(%) 



Statistical eiTor in calibration ratio 0.60 

^-dependence of calibration ratio 1.10 

Uncertainty in B03 Bi 1.10 

Uncertainty in QUaD Bg 0.75 

Pointing uncertainty 1 .62 

Internal consistency (jackknifes) 1.20 

B03 calibration error 2.00 

Total uncertainty 3.38 



To propagate the errors in the B03 and QUaD beam transfer functions onto our calibration, we repeat the analysis but with the 
beam functions shifted by their quoted errors. Doing this for each of the B03 and QUaD beams, we take the resulting shifts in 
the calibration numbers as the error due to uncertainty in the beams. We find a 1.1% shift due to the uncertainty in the B03 beam 
function and a 0.75% shift due to the uncertainty in QUaD's beam. 

A further source of error is the relative pointing uncertainty between the QUaD and B03 maps. There is a clear pointing offset 
seen between the QUaD and B03 maps and so we have shifted the B03 maps before performing the calibration analysis. We find 
the appropriate shift (which we model as a simple shift in R.A. and Dec) by fitting for it in map-space. To quantify the error, 
we repeat the analysis with the B03 maps shifted (away from the best fit) according to the errors returned from our map-based 
fit. Applying the shift in a number of different directions (the eight compass points), we find the maximum shift in the resulting 
calibration factor is 1.6%. We take this number as our error due to the B03/QUaD relative pointing uncertainty. 

Finally, QUaD also inherits the stated uncertainty in the B03 calibration which is 2% (Masi et al. 2006). We add this and each 
of the errors derived above in quadrature to arrive at our final calibration uncertainty of 3.4%. Our absolute calibration error 
budget is summarized in Table |7] 
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